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Abstract. We have studied the performance of global fitting 
of low-resolution X-ray spectra in retrieving intrinsic source 
parameters, with emphasis on the coronal metallicity. The 
study has been conducted by fitting large numbers of simulated 
spectra with known characteristics, and studying the distribu- 
tion of best-fit parameters. We have studied the behavior of the 
LECS detector on board the SAX satellite and the SIS detector 
on board the ASCA satellite. The fitted source spectra have ei- 
ther two discrete temperature components or a power-law tem- 
perature distribution, with metallicity variations modeled by a 
single global abundance parameter. The model used for the fit- 
ting has always been a two-temperature one, with global vary- 
ing abundance, to explore the influence of the a priori ignorance 
of the actual temperature stratification in the source being ob- 
served. 

The simulations performed explore the influence of vary- 
ing statistics in the observed spectrum (spanning a realistic 
range of values) as well as the effect of varying the intrin- 
sic source metallicity, with values in the range 0.15-1.0 times 
the solar value. We find that the source metallicity can be re- 
trieved within few tens of percent from ASCA/SIS spectra 
of typical signal to noise ratio, and within few percent from 
SAX/LECS spectra at the same signal to noise ratio. However 
relatively small uncertainties in the detector calibrations and in 
the plasma emission codes are likely to potentially cause large 
systematic off-sets in the value of the best-fit parameters. Simi- 
lar systematic off-sets may derive from assuming too simplistic 
a temperature distribution for the source plasma. 

In addition we have re-analyzed the ASCA/SIS spectra of 
the active giants /? Cet and Capella with the same set of as- 
sumptions used in the simulations, showing how the best-fit 
metallicity in these two real cases depends on the details of the 
fitting process, and in particular on the chosen energy range. 

Key words: Stars: abundances; stars: late-type; X-rays: stars; 
methods: data analysis 
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1. Introduction 

CuiTently available soft X-ray spectrographs are of the non- 
dispersive type (CCDs and GSPCs), and the "analysis" of 
a spectrum, as commonly performed, consists essentially in 
searching, given a detector response matrix and a plasma emis- 
sion code, the space of possible "models" (i.e. distributions of 
temperatures and relative normalizations, as well as elemental 
abundances), for the "best-fit" one. This is usually done within 
a set of a priori constraints on the model, in particular fixing 
the number of discrete temperature components allowed. The 
usual measure of "goodness of fit" used is the x^, and the fit- 
ting process thus consists in a search, through space, for the 
model yielding the minimum value. 

One notable result of X-ray coronal astronomy coming 
from ASCA is that most (if not all) coronal spectra observed 
so far with the SIS CCD detector yield, when fit with currently 
available plasma emission codes, best-fit parameters incompat- 
ible with solar (or near-solar) abundance plasmas. At the same 
time they appear to be "better fit" with elemental abundances 
typically a fraction of the solar abundance, with values as low 
as 0. 1 times solar not uncommon. This finding points toward 
a widespread lower-than-solar metal abundance in stellar coro- 
nae, although departures from the stellar photospheric abun- 
dances, which would be of greater importance, have not been 
usually investigated. 

The limited spectral resolution of non-dispersive X-ray de- 
tectors makes it impossible to determine reliably the flux of 
individual lines or line complexes, so that most studies of coro- 
nal abundances up to now have been conducted by doing global 
fits to the whole spectrum, using the resulting best-fit parame- 
ters (including the abundances) as the most probable source 
parameters. In the presence of complex spectral models, how- 
ever, it is difficult to assess from first principles what the actual 
uncertainties related to the fitting process are, and how are the 
various fit parameters (possibly) correlated with each other. 
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To help assess the significance of the abundance estimates 
derived from low-resolution X-ray spectra, as well as to help 
compare the relative performance of the SIS and LECS instru- 
ments for this type of studies, we have performed an exten- 
sive set of simulations, to study the sensitivity of the two in- 
struments to variations in the metallicity in the source spec- 
trum, and to assess the relative uncertainties in the derived 
source parameters. We have first assumed that the source spec- 
trum is intrinsically emitted by two distinct isothermal com- 
ponents, i.e. a so-called "two-temperature" model. This model 
is bound to be a simplistic representation of a more compli- 
cated reality, with the coronal plasma of real stars certainly 
having a much more complex distribution of temperatures, but 
the two-temperature model is nevertheless still the work-horse 
of present day non-dispersive X-ray coronal spectroscopy. Pre- 
vious work has shown that the spectrum of multi-temperature 
plasmas confined in coronal loops and yielding a moderate 
number of counts (as the ones discussed here) may be rea- 
sonably well fit by two-temperature models (CiaraveUa et al. 
1997). This evidence justifies, for the purpose of the present 
paper, the usage of two-temperature models as "representative" 
of a more complex reahty. The two-temperature source spectra 
have been fitted with two-temperature models. 

As a check for the dependence of the fitting process on the 
correctness of the assumed source spectrum we have performed 
an additional set of simulations assuming that the source has a 
power-law temperature distribution, with spectral index a = 
1.5 and maximum temperature of 3.0 keV. This type of power- 
law distribution is reasonably representative, for example, of 
the emission measure distribution of the solar corona in the 
temperature range accessible to soft X-ray detectors (Raymond 
& Doyle 1981). These power-law spectra were fit with the same 
type of two-temperature models as the two-temperature source 
spectra, thus simulating our ignorance about the "true" intrinsic 
source temperature distribution. In the present work the plasma 
emission model used has been the MEKAL code (Mewe et al. 
1995), as implemented in the XSPEC X-ray spectral analysis 
package (version 9.0) 

2. Characteristics of the instruments 

The LECS instrument (Parmar et al. 1997) is a drift-less gas 
scintillation proportional counter with a thin entrance window, 
providing continuous energy coverage in the spectral range 
from 0.1 to lOkeV. Its resolution is energy-dependent (vary- 
ing as E~^-^), and becomes comparable to the resolution of 
the CCD-based SIS detector at ~ 0.5 keV. 

The SIS instrument is a CCD-based non-dispersive X-ray 
spectrometer, providing energy coverage in the band from 0.5 
to lO.OkeV, with an approximately constant energy resolution, 
which is effectively decreasing with time because of the radia- 
tion damage of the CCD chip (Dotani et al. 1996). At the time 
of the launch the effective resolution was ~ 100 eV at 1 keV, 
with the precise value depending on the details of the instru- 
mental mode being used. However, the number of spectral res- 
olution elements for the SIS and for the LECS is similar, due 



to the larger passband of the latter. The spectral resolution of 
the LECS and of the SIS (at launch) are compared in Fig. 1. 
We have used a SIS response appropriate for an observation 
made in early 1995, i.e. two years into operations. This implies 
a detector performance intermediate between the best possible 
performance (at launch) and present-day performance. 

The effective area of the SIS instrument is, thanks to its 
bigger mirror, larger than the effective area of the LECS, thus 
giving it access to fainter sources. The wider passband of the 
LECS partially compensates, for soft source spectra as the one 
discussed here, the smaller effective area of the mirrors, so that 
a given source with the spectral parameters used here yields a 
comparable count rate both in the SIS and in the LECS. The ef- 
fective areas of the SIS and of the LECS detector are compared 
in Fig. 2. 




0.1 1 10 

Energy (keV) 



Fig. 1. The full width at half maximum (FWHM) spectral resolution 
of the SAX/LECS (continuous line) and ASCA/SIS (dashed hne) in- 
struments. The SIS resolution is the one at launch time. 



3. The simulations 

In the first set of simulations a two-temperature source spec- 
trum was used, with temperatures of 0.5 and 2.0 keV, and iden- 
tical emission measure. Real coronal sources, when modeled 
with a two-temperature spectrum, display a wide range of pa- 
rameters, but the values assumed here are reasonably repre- 
sentative of "typical" active stars, as determined, for example, 
by fitting SIS data. Since the "true" coronal source spectra are 
likely to be produced by a plasma with a more complex tem- 
perature and density distribution, as inferred from spatially- 
resolved observations of the solar corona, a second set of simu- 
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Fig. 2. The effective areas of the SAX/LECS (continuous line) and 
ASCA/SIS (dashed hne) instruments. 



lations was performed, using a power-law emission measure 
distribution with slope a = 1.5 and maximum temperature 
Tmax = 3 keV for the source spectrum. These parameters are 
comparable to those found by Preibisch (1997) when fitting 
ROSAT/PSPC spectra of field solar-type stars. 

The XSPEC package has been used both to produce the 
simulated spectra and to fit them. 300 realizations of the source 
spectrum with fixed parameters and with a fixed total number 
of counts, including statistical errors determined by Poisson 
statistics were generated. The simulated spectra were then re- 
binned, so to have at least 20 counts in each resulting new (vari- 
able size) energy bin, to ensure that the statistics is appUca- 
ble (Eadie et al. 1971), and each spectrum was fitted, always 
with a two-temperature mekal model, retrieving the best-fit 
spectral parameters. 

The simulations assuming a two-temperature source spec- 
trum were run with 2500, 10000 and 40000 source counts, 
spanning the typical source counts for actual SIS and LECS 
spectra. The simulations assuming an intrinsic power-law spec- 
trum were all performed with 10000 source counts. Simula- 
tions can supply at most a lower limit to the uncertainty in the 
derived spectral parameters, because of at least two optimistic 
assumptions which are only approximations to the reality. First, 
we are assuming that the mekal plasma emission code mod- 
els describe perfectly the emission from a real plasma, neglect- 
ing, among other things, all of the uncertainties in the atomic 
physics. Second, we are assuming that X-ray detectors are per- 
fectly calibrated, neglecting (possible) systematic uncertainties 
present in the whole process. 



Table 1. The numerical experiments performed for the present paper. 
Two-temperatures source spectra were generated with temperatures of 
0.5 and 2.0 keV, and with an emission measure ratio of 1. Power-law 
spectra have a slope of —1.5 and a maximum temperature of 3.0 keV. 



\_oae 


counts 


z 


Detector 




2-T source 


spectra 


A2TL0.15 


2500 


0.15 


ASCA/SIS 


A2TM0.15 


10000 


0.15 


ASCA/SIS 


A2TH0.15 


40000 


0.15 


ASCA/SIS 


A2TL0.50 


2500 


0.50 


ASCA/SIS 


A2TM0.50 


10000 


0.50 


ASCA/SIS 


A2TH0.50 


40000 


0.50 


ASCA/SIS 


A2TL1.00 


2500 


1.00 


ASCA/SIS 


A2TM1.00 


10000 


1.00 


ASCA/SIS 


A2TH1.00 


40000 


1. 00 


ASCA/SIS 


S2TL0.I5 


2500 


0.15 


SAX/LECS 


S2TM0.15 


10000 


0.15 


SAX/LECS 


S2TH0.15 


40000 


0.15 


SAX/LECS 


S2TL0.50 


2500 


0.50 


SAX/LECS 


S2TM0.50 


10000 


0.50 


SAX/LECS 


S2TH0.50 


40000 


0.50 


SAX/LECS 


S2TL1.00 


2500 


1.00 


SAX/LECS 


S2TM1.00 


10000 


1.00 


SAX/LECS 


S2TH1.00 


40000 


1.00 


SAX/LECS 


power-law EM source spectra 


APLM0.15 


10000 


0.15 


ASCA/SIS 


APLM0.50 


10000 


0.50 


ASCA/SIS 


APLMl.OO 


10000 


1.00 


ASCA/SIS 


SPLM0.15 


10000 


0.15 


SAX/LECS 


SPLM0.50 


10000 


0.50 


SAX/LECS 


SPLMl.OO 


10000 


1.00 


SAX/LECS 



Additionally, we are assuming that the two-temperature 
model or the power-law emission measure distribution are ad- 
equate representations of the source spectrum. On the basis of 
our knowledge of the solar corona, as well as on the basis of 
analyses of EUV spectra of real coronal sources, this is likely 
to be, at best, an approximate parameterization of the coronae 
of real stars. 

Abundance variations in the source spectrum have been 
modeled through a single parameter, the global metalUcity Z, 
thus assuming no variations in the abundance ratios with re- 
spect to the solar values. In fact, several of the analyses pub- 
Ushed on SIS spectra of coronal sources result in best- fit spectra 
with significant variations in the elemental abundance ratios. A 
simulation-based analysis for the case of individually varying 
elemental abundances will be the subject of a future paper. 

Each numerical experiment is indicated in the following, 
with a code, and the correspondence between the codes and the 
experiment parameters is given in Table 1 . 
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4. Simulation results 



Each simulation yields a set of best-fit spectral parameters for 
all the realizations of the sample spectrum. The distribution of 
the best-fit values is indicative of the uncertainties that can be 
expected, in the absence of systematic effects, when fitting real 
spectra with similar statistics. 

Some representative cases for the simulations assuming 
two-temperature source spectra are shown in Figs. 3 and 4, in 
the form of scatter diagrams, plotting best-fit spectral parame- 
ters against each other. For each simulation we have computed 
the median, as well as the upper and lower 68% and 90% quan- 
tiles, corresponding to "1 cr" and "2ct" levels for each best-fit 
parameter. The median and quantiles are shown in numerical 
form in Table 2, and in graphical form in Figs. 5 and 6. 
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Fig. 3. Two-temperature source spectra fit with two-temperature mod- 
els: distribution of the best-fit parameters for a set of 300 realiza- 
tions, with coronal abundance 0.5 times solar as observed with the 
SAX/LECS detector, with 10000 counts accumulated per spectrum 
(experiment S2TM0.50). 



4.1. Results for the SAX/LECS two-temperature simulations 

Our simulations show that the combination of resolution and 

spectral coverage offered by the LECS detector allows retrieval 
of the spectral parameters through a fit, for two-temperature 
spectra, with good accuracy, already for spectra with ~ 10 000 
counts The 2-a accuracy is ~ 5% for the source temperatures, 
and ~ 15% for the coronal metallicity. LECS spectra with 
~ 40 000 counts allow 2-a accuracies of ~ 2% on the source 
temperatures, and ~ 7% on the coronal metallicity. 

The best-fit parameters show some correlation (Fig. 3), i.e. 
they appear to be mildly dependent on each other in statistical 
terms, but the range of variation of the metallicity vs. the emis- 
sion measure ratio or the temperatures is relatively small. The 
scatter plots for the rest of the simulations are very similar to 
the ones in Fig. 3, with the same qualitative features, and a to- 



Fig. 4. Two-temperature source spectra fit with two-temperature mod- 
els: distribution of the best-fit parameters for a set of 300 realiza- 
tions, with coronal abundance 0.5 times solar as observed with the 
ASCA/SIS detector, with 10000 counts accumulated per spectrum 
(experiment A2TM0.50). 



tal scatter decreasing as expected with increasing total number 
of counts, as summarized in Fig. 5. At the same time, as shown 
in Fig. 6, the relative scatter is essentially independent on the 
intrinsic source metalUcity. 

4.2. Results for the ASCA/SIS two-temperature simulations 

The results for the simulations relative to the SIS detector dif- 
fer from the ones relative to the LECS, in the case of two- 
temperature source spectra fit with two -temperature models, in 
one important way: the scatter of the best-fit coronal metallicity 
value is (for the same total number of counts in the spectrum) 
much larger, as it is the scatter for the best-fit emission-measure 
ratio. The two are also strongly correlated with each other. 

The larger observed scatter implies that while the temper- 
atures can be measured in a ~ 10 000-counts SIS spectrum 
with a 2-(7 accuracy of ~ 5%, the coronal metalhcity can be 
measured, in the same spectrum, with a 2-cr accuracy of only 
~ 40%, and an accuracy on the emission-measure ratio of 
~ 25%. Even in the higher statistics spectra (~ 40 000 counts), 
the 2-cr accuracy for the source metallicity is still only ~ 20%. 
This, again, in the optimistic assumption of an intrinsic two- 
temperature source spectrum. 

Again, the scatter plots not shown have the same qualitative 
features of the ones shown, with the total scatter decreasing, as 
expected, with increasing total number of counts. 

4.3. Results for the SAX/LECS power-law simulations 

The confidence ranges for the best-fit parameters obtained by 
fitting a power-law source spectrum with a two-temperature 
model in the case of the LECS detector are shown in Fig. 7, 
together with the results relative to the same type of simula- 
tions done for the SIS. In this case the confidence ranges are 
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Table 2. Observed 1- and 2-a ranges for the best-fit parameters for all the numerical experiments discussed in the present paper. The entries in 
the table are in percent of the median value of the relevant best-fit parameter. 



Code Ti T2 EM1/EM2 Z 





1 

i- 


■G 


L- 


(J 


1 ^ 




L- 


•G 


1 

i- 


G 


Z-i 


tj 


1 

1- 


G 




G 


A2TL0.15 


+7 


— 7 


+13 


— 15 


+13 — 


6 


+29 


— 12 


+74 


— 29 


+175 


— 40 


+71 


— 44 


+114 


— 65 


A2TM0.15 


+3 


_ 3 


+5 


- 5 


+2 - 


4 


+5 


_ 7 


+27 


- 13 


+49 


— 23 


+24 


- 25 


+46 


- 34 


A2TH0 1 5 


-1-1 


_ 1 




_ 3 




3 




— 4 


+12 


_ 5 


4-21 


— 10 


+7 


— 13 


+14 


— 21 


A2TL0.50 


-1-4 


— 4 


+7 


— 8 


+11 — 


7 


+18 


— 11 


+35 


— 23 


+70 


— 37 


+59 


— 32 


+131 


— 46 


A2TM0.50 


+2 


- 2 


+3 


- 4 


+3 - 


4 


+6 


- 7 


+15 


- 11 


+29 


- 18 


+25 


- 15 


+45 


- 29 


A2TH0.50 


+ 1 


_ 1 


+1 


- 2 


+0 - 


3 


+1 


- 4 


+10 


- 2 


+ 14 


- 7 


+6 


- 11 


+14 


- 16 


A2TL1.00 




_ 2 


+4 


— 4 


+17 — 


7 


+26 


— 14 


+36 


— 24 


+86 


— 36 


+105 


— 32 


+214 


— 58 


A2TM1.00 


+ 1 


- 2 


+3 


- 3 


+3 - 


5 


+6 


- 7 


+17 


- 9 


+27 


- 16 


+28 


- 21 


+52 


- 31 


A2TH1.00 


+1 


- 1 


+1 


- 1 


+0 - 


3 


+2 


-4 


+9 


-3 


+14 


-8 


+10 


- 11 


+21 


- 18 


S2TL0.15 


+11 


— 11 


+16 


— 20 


+13 — 


5 


+20 


— 11 


+37 


— 15 


+58 


— 28 


+21 


— 19 


+58 


— 28 


S2TM0.15 


+5 


- 5 


+9 


- 9 


+4 - 


4 


+7 


— 6 


+18 


- 6 


+25 


— 13 


+10 


- 9 


+19 


- 15 


S2TH0.15 


+3 


- 2 


+5 


- 4 


+0 - 


3 


+1 


- 4 


+6 


- 4 


+ 12 


- 7 


+5 


- 5 


+9 


- 8 


S2TL0.50 


+7 


- 5 


+10 


- 12 


+11 - 


5 


+19 


- 10 


+26 


- 12 


+41 


- 19 


+21 


- 17 


+42 


- 24 


S2TM0.50 


+3 


- 2 


+5 


-5 


+4 - 


3 


+6 


-5 


+13 


-4 


+22 


-9 


+8 


- 7 


+16 


- 12 


S2TH0.50 


+1 


- 1 


+2 


-3 


+0 - 
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+2 


-3 


+7 


- 2 


+10 


-4 


+4 


-4 


+6 


-7 


S2TL1.00 


+5 


- 5 


+8 


- 9 


+10- 


6 


+18 


- 10 


+20 


- 12 


+35 


- 18 


+26 


- 17 


+49 


- 25 


S2TM1.00 


+2 


- 2 


+4 


- 4 


+3 - 


3 


+6 


-5 


+12 


- 2 


+17 


- 7 


+10 


- 7 


+17 


- 14 


S2TH1.00 


+1 


- 1 


+9 


-7 


+3 - 


3 


+5 


-5 


+6 


-0 


+9 


-3 


+3 


-4 


+7 


-7 


APLM0.15 


+3 


-3 


+5 


- 5 


+4 - 


4 


+6 


-6 


+9 


- 8 


+15 


- 13 


+16 


- 14 


+29 


- 21 


APLM0.50 


+4 


- 4 


+6 


- 7 


+4 - 


3 


+6 


-6 


+15 


- 10 


+23 


- 14 


+14 


- 15 


+29 


- 24 


APLMl.OO 


+10 


-8 


+37 


- 11 


+4 - 


4 


+26 


-6 


+36 


- 22 


+403 


-30 


+22 


-24 


+33 


-60 


SPLM0.15 


+6 


- 4 


+9 


- 7 


+3 - 


3 


+5 


- 5 


+18 


- 13 


+33 


- 21 


+8 


-9 


+ 14 


- 15 


SPLM0.50 


+6 


- 7 


+13 


- 11 


+4 - 


3 


+9 


-5 


+26 


- 17 


+55 


-28 


+11 


- 10 


+20 


- 17 


SPLMl.OO 


+21 


- 16 


+30 


-29 


+19- 


8 


+32 


- 12 


+130 


-43 


+234 


-59 


+33 


- 19 


+57 


-28 



not computed relative to the "true" source parameters (as the 
source model is different from the model used to fit the data, 
and thus there are no "true" two temperatures for the source) 
but rather with respect to the median value found in the simu- 
lated set. The only parameter for which the confidence regions 
are computed relative to the source value is the metal abun- 
dance. Inspection of Fig. 7 and of Table 2 shows that, if the 
intrinsic source metallicity is really as low as 0.15 times so- 
lar, the confidence regions derived by using a two-temperature 
model are exceedingly broad. We will further discuss only the 
results relative to the Z = 1.0 and Z = 0.5 cases. 



A comparison of Fig. 6 with Fig. 7 shows that, for the 
LECS detector, the size of the confidence regions for the best- 
fit temperatures has increased slightly in comparison with the 
case of two-temperature source spectra, as it has for the best- fit 
emission-measure ratio. The best-fit metallicity is also show- 
ing slightly larger confidence regions, in addition to a signifi- 
cant offset from its true value, with the two-temperature model 
retrieving lower abundances than the true value. 



4.4. Results for the ASCA/SIS power-law simulations 

When the SIS is used to fit power-law source spectra with a 
two-temperature model, the on average does not converge 
to a satisfactory value. For the cases with Z = 0.5 and Z = 1.0 
the median values of the fit are 1.17 and 1.27, respectively, 
which on account of the large number of degrees of freedom of 
the fits presented here (~ 200) corresponds to fits which are 
formally unacceptable with a probability (for the Z = 1 case) 
of ~ 0.5%. This has an influence on the size of the confidence 
regions, and makes it appear the fitting process to be "more pre- 
cise" than it actually is. As a result, the confidence regions for 
this case cannot be directly compared with the previous cases, 
in all of which the best fit distribution did not show any 
excess above the expected distribution. 

The confidence regions for the SIS two-temperature fits 
of power-law source spectra are still significantly larger than 
the LECS ones, although they appear (because of the above- 
mentioned effect) to be smaller than the regions for the two- 
temperature fitting of two-temperature source spectra (again 
but for the case Z = 0.15, which shows an anomalous scatter). 
However, the median best-fit metallicity appears to be close to 
the intrinsic metaUicity of the original power-law model. 
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Fig. 5. Two-temperature source spectra fit with two-temperature mod- 
els: the effect of increasing photon statistics on the 68% and 90% error 
ranges for the best-fit parameters. Each panel shows the error ranges 
for a given parameter, with the number of counts in the total spectrum 
indicated by the bottom label. The top label indicates the instrument 
("S" for SAX/LECS and "A" for ASCA/SIS). The two top panels show 
the error range for the cool (top left) and hot (top right) temperatures, 
while the bottom panels show the range for the ratio between the two 
emission measures (bottom left) and for the coronal metallicity (bot- 
tom right). 

5. Comparison between the ASCA/SIS and SAX/LECS de- 
tectors 

The SIS and LEGS detectors differ in their trade-offs of resolu- 
tion and spectral coverage, and thus behave differently in their 
capability to retrieve source parameters, as well as in their sen- 
sitivity to discrepancies between the intrinsic source spectral 
model and the model being used to fit the data. 

5.7. Two-temperature intrinsic source spectra 

Given the same total number of counts in the spectrum, the SIS 
detector does a slightly better job in determining the source 
temperatures, but it performs worse than the LEGS in terms 
of retrieving the coronal metallicity of the source under inves- 
tigation (Figs. 5 and 6). A solar-metallicity, ~ 10 000-counts 
LEGS spectrum allows the coronal metalUcity to be determined 
with a 2-cr accuracy of ±15%, comparable to the accuracy ob- 
tained with a ~ 40 000-counts SIS spectrum. 

The reason for this large difference Ues in the diagnostic 
which is being exploited by the fitting process to derive the 
abundance. The metallicity is essentially determined by balanc- 
ing the spectrum in the region between 0.7 and 1.2keV, where 
Une emission (mostly by FeL lines) dominates and where the 



Fig. 6. Two-temperature source spectra fit with two-temperature mod- 
els: the effect of different source metallicity on the 68% and 90% er- 
ror ranges for the best-fit parameters, illustrated using the simulations 
with 10 000 counts in the source spectrum. Each panel shows the er- 
ror ranges for a given parameter, with the source metallicity indicated 
by the bottom label. The top label indicates the instrument ("S" for 
SAX/LECS and "A" for ASCA/SIS). The panels have the same struc- 
ture as in Fig. 5. 



statistics are best, both in the SIS and in the LEGS, with the 
available continuum emission. The region around 1 keV is thus 
acting as a sort of "pivot" for the fitting process. In the case of 
the SIS, Uttle continuum is usually available, mostly in the (lim- 
ited) region between 0.5 and 0.7 keV, in which however the de- 
tector's resolution is lowest, and thus the continuum contribu- 
tion may not be easily disentangled from the nearby line emis- 
sion. The continuum in the hard tail of the spectrum, typically 
the small regions between the major K-complexes of heavier 
elements has, for most coronal sources, low intrinsic source 
flux, and thus low signal to noise ratio. 

For LEGS spectra the region below 0.5 keV supplies a very 
well constrained determination of the continuum, and makes 
the fitting process much more robust. To confirm the impor- 
tance of the region below 0.5 keV we have performed an ad- 
ditional set of simulations for the LEGS, fitting only the spec- 
trum above 0.5 keV, i.e. using only the information available to 
the SIS detector. As expected, the performance for the LEGS 
becomes comparable to the performance of the SIS (actually 
sUghtly worse, not surprisingly, given its lower spectral res- 
olution). Fig. 9 shows the effect of moving the low-energy 
cut from 0.1 to 0.5 keV for LEGS spectra on the confidence 
regions, for the sample case of solar-metalUcity spectra with 
10000 counts. 
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Fig. 7. Power-law emission measure source spectra fit with 
two-temperature models: the effect of different source metallicity on 
the 68% and 90% error ranges for the best-fit parameters. Each panel 
shows the error ranges for a given parameter, with the source metal- 
licity indicated by the bottom label. The top label indicates the instru- 
ment ("S" for SAX/LECS and "A" for ASCA/SIS). The panels have 
the same structure as in Fig. 5. 
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Fig. 8. The two-temperature model source spectrum, for solar abun- 
dance, used in the present work. Notice the almost line-free continuum 
below 0.5 keV, which can be exploited by the SAX/LECS for metal- 
licity determinations and which is not available to the ASCA/SIS. 
The two dashed lines are the two individual temperature components, 
while the continuos line is the total spectrum. 



Fig. 9. Two-temperature source spectra fit with two-temperature mod- 
els: the effect of moving the low-energy cut for SAX/LECS spec- 
tra from 0.1 to 0.5 keV is illustrated for the sample case of so- 
lar-metallicity spectra with 10 000 counts. The panels have the same 
structure as in Fig. 5. 

5.2. Power-law intrinsic source spectra 

When power-law source spectra are fit with two-temperature 
models, the two detectors discussed here behave differently. 
While the LECS produces formally acceptable fits, but with 
lower best-fit metalUcities than the intrinsic source value, the 
SIS detector produce formally unacceptable fits, while still 
converging on the right source metallicity. 

This behavior can be explained as follows: when perform- 
ing simulations with the same total number of counts for the 
SIS and the LECS, the smaller energy coverage of the SIS re- 
sults in a concentration of counts in the ~ 1 keV region, yield- 
ing spectra with higher statistics than in the case of the LECS, 
in which the counts are spread in a larger spectral region with 
two distinct peaks (at ~ 0.2 and ~ 1 keV), so that the result- 
ing statistics (particularly in the peak channels) are correspond- 
ingly lower. This explains why the SIS spectra in this case re- 
sult, with the same number of total counts, in worse values 
than the LECS ones. 

The convergence of LECS spectra on a lower than intrin- 
sic metalUcity value is due to a different effect: LECS fits are 
driven mainly by the balance between the two high-statistics 
regions of the spectrum, i.e. the continuum below ~ 0.4 keV 
and the Une-rich region around ~ 1 keV. When attempting to 
fit power-law spectra with two-temperature models the discrep- 
ancy between the source spectrum and the fitting model is such 
that the fit is lead to a lower abundance model by the need to 
balance the Unes and the continuum. This does not happen in 
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the case of the SIS, again because of the lack of high- statistics 
continuum regions. 

The reality of the above effect is well shown by a further 
set of simulations, in which the same LEGS power-law source 
spectra are fit only using the spectral region above 0.5 keV. In 
this case the fit is "pivoting" around the high-statistics 1 keV 
region, and is being driven to higher metallicities (this time 
higher than the intrinsic source metallicity) by the residuals 
in the high-energy tail. Fig. 10 shows how the confidence re- 
gions change when the soft cut in the LEGS spectra is moved 
from 0. 1 to 0.5 keV, indicating that when the fitting model is not 
strictly adhering to the true source thermal distribution the best- 
fit metallicity can be strongly influenced by the spectral region 
being used. Conversely, a best-fit metallicity which changes 
with the spectral region being fit is simptomatic of too sim- 
pUstic a model. The direction and the size of the effect in real 
cases will depend on the (a priori unknown) intrinsic source 
temperature distribution, so that the present case only has an 
illustrative value. 
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Fig. 10. Power-law emission measure source spectra fit with 
two-temperature models: the effect of moving the low-energy cut for 
SAX/LECS spectra from 0. 1 to 0.5 keV is illustrated for the sample 
case of solar-metallicity spectra with 10 000 counts. 



6. Two test cases: the observed ASCA/SIS spectra of (3 Ceti 
and of Capella 

We have analyzed the SIS spectra of the active giant /3 Get and 
of the active binary GapeUa obtained during the Performance 
Verification (PV) phase. The data were retrieved from the pub- 



lic archive, and analyzed according to the procedures suggested 
for ASGA data. 

The data for /3 Get consist of two segments, of ~ 8 and ~ 
3.4 ks each, taken in two different instrumental modes (faint 
and BRIGHT mode). They have thus been processed separately, 
but have been fit simultaneously, as discussed in Maggio et al. 
(1997). Although the /3 Get SIS spectinm is better fit with a 
model in which the individual abundances are left free to vary 
(as is the GapeUa spectrum), for coherence with the simula- 
tions, they have been fit here with global varying abundances. 
Again, we have used the MEKAL emission code as imple- 
mented in XSPEG V.9.0; the results therefore are not directly 
comparable with the ones of Drake et al. (1994), obtained with 
the MEKA emission code (an earlier version of the MEKAL 
code, see Mewe et al. 1995 for a description of the changes). 
Eventual differences in the version of the ASGA pipeline pro- 
cessing may also influence the results. The two resulting spec- 
tra have been rebinned at a minimum of 20cts per bin. The 
summed number of counts from both spectra was ~ 10 000. 

The f3 Get SIS spectrum was first fitted across a broad band, 
from 0.5 keV up to 4.0 keV. In this energy range the best-fit 
metallicity is ~ 0.9 times solar. However, the exclusion of 
channels at the soft end of the spectrum causes the best-fit 
metallicity to change, making it decrease to 0.33 solar if chan- 
nels with energies from 0.5 to 0.6 keV are excluded, and as 
low as 0.22 if channels up to 0.7 keV are excluded. At the 
same time, the ratio between the emission measures of the cool 
and of the hot component increases, showing the same type 
of behavior as the simulations for the case of intrinsic two- 
temperature source spectra, where lower best-fit metallicities 
are correlated with higher emission-measure ratios. 

The same type of analysis was done for the Gapella spec- 
trum, finding a best-fit abundance for the broad-band (0.5- 
5.0 keV) SIS Gapella PV phase spectrum of 0.54 times solar. 
As for (3 Get, the exclusion of spectral channels below 0.6 keV 
lowers the best-fit abundance to 0.24 times solar, and exclu- 
sion of the channels below 0.7 keV lowers it even more, down 
to 0.18 times solar. As for /? Get these results are not directiy 
comparable to the ones of Drake (1996). 

The does not converge, neither for GapeUa, nor for /3 Get 
to formally acceptable values. The physical meaning of best-fit 
parameters is thus clearly questionable, and they should not 
therefore be taken at face value. The strong and systematic de- 
pendence of the best-fit metallicity on the spectral region in the 
Gapella and /? Get spectra indicates the presence of some sys- 
tematic error either in the spectral model used (i.e. the emission 
code itself, the source temperature structure, the assumption of 
solar metallicity ratio or a combination of these) or in the de- 
tector's caUbration (or both). 

The continuum region used by the SIS fits (the region be- 
tween ~ 0.5 and ~ 0.75 keV) is almost (but not entirely) line- 
free, but, at the resolution of the SIS, and given the caUbration 
uncertainty in the region (Dotani et al. 1996) it contains a sig- 
nificant number of source counts "spilling over" from the Fe-L 
complex. The statistics in this region are, for coronal sources, 
invariably higher than for the hard tail where line complex from 
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the K-transitions of the heavier elements are visible, and will 
accordingly influence the metallicity determination. 

The increase of the best-fit metalhcity with inclusion of a 
softer spectral region can be explained by assuming that the 
models being used for the fit (either because of problems in 
the emissivity models or because of cahbration problems, or a 
mixture of both) have a lack of flux in the region around 1 keV 
(which is dominated by the cooler temperature component). 
This could happen, for example, if a number of lines which 
exist in real spectra were missing from the models. In this case 
the fit process would converge by balancing the apparent flux 
deficit around the line-rich 1 keV region. The strategy "chosen" 
by the fit process will be different depending on the spectral 
region being fit. If the softer channels are excluded, the fit pro- 
cess can freely increase the cool emission measure to drive up 
the 1 keV region, while at the same time decreasing the metal- 
licity to avoid discrepancy in the hard tail, dominated by the 
hot component. If the softer channels are included, however, 
the high emission measure of the cool component drives up the 
pseudo-continuum around 0.5 keV, and thus the fit process has 
to increase the metalhcity (to add flux around 1 keV without 
adding more continuum at 0.5 keV) and correspondingly de- 
crease the emission measure of the soft component. The effect 
is well visible in the SIS spectrum of (3 Cet shown in Fig. 1 1 . 

We have also fit the same spectra ignoring the region be- 
tween 0.75 and 1.1 keV. The result reinforces the above con- 
clusions: the best-fit metallicity for the /? Cet spectrum rises to 
1.1 times solar, the same value which is found for the Capella 
spectrum applying the same procedure. The lower metaUicity 
obtained by fitting the whole spectrum could be explained with 
a lack of flux (again due to calibration effects and/or to lack of 
Unes in the emission models) around the 1 keV region. Another 
possibility (which is not mutually exclusive with the first one) 
is that the two-temperature model being used to fit the data is 
not a good representation of the real temperature distribution 
of the emitting plasma. 

The possible lack of lines in current plasma emission codes 
in the region around IkeV has for example been discussed 
by Jordan (1996), who makes explicit reference to the MEKA 
model (from which the MEKAL model used here has been de- 
rived by improving the Fe ionization balance and by adding a 
number of lines), concluding that "in view of the large number 
of weak transitions not yet included in the emissivity codes it 
seems premature to attribute discrepancies between the obser- 
vations and the predictions of the codes to non-photo spheric 
abundances" . 

The dependence of derived abundances on the precise spec- 
tral interval being fit does not disappear if a model with all 
the abundances free to vary is used. Rather, it becomes more 
complicated because of the interplay between the various abun- 
dances, and because introducing a cut-off in the soft spectral 
region will make some specific abundances impossible to de- 
termine (notably nitrogen and oxygen). However, the resulting 
metaUicities stiU strongly depend on the spectral interval being 
considered. 
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Fig. 11. The ASCA/SIS 13 Cet spectra, together with the best-fit model 
determined only considering the spectral range above 0.7 keV. Note 
that this model significantly over-shoots, as described in the text, the 
data in the softer region 0.5-0.7 keV. 

7. Discussion 

The set of simulations presented here allow us to get some 
insight in the complex process of determining coronal abun- 
dances by global fitting of X-ray spectra, as well as to deter- 
mine the expected range of uncertainty for the best-fit param- 
eters as a function of spectrum statistics and of the intrinsic 
source metallicity. In addition, by giving an estimate of the 
range of variation in best-fit parameters that can be induced 
by a certain level of statistical error, it allows to estimate what 
can be the (approximate) effect of a certain level of systematic 
error present, for example, in the detector calibration or in the 
emissivity model. 

The LECS performs better in constraning the intrinsic 
source metaUicity than the SIS. This is due to the larger spectral 
range accessible, and in particular to the avaUabiUty of a high 
signal-to-noise continuum region below 0.5 keV, with few and 
relatively weak lines, which appears to be more important, for 
spectroscopy of coronal sources, than the enhanced resolution 
of the SIS. In fact, as reported in Sect. 5.1, when we remove 
the spectral region below 0.5 keV from the LECS spectra, thus 
making its passband equal to the one of the SIS, the indeter- 
mination on the metallicity derived from spectra of the two 
instruments is similar, as it is the strong correlation between 
metallicity and emission measure ratio. 

Some results from the simulations are remarkably differ- 
ent from what could have been expected. For example the cool 
component of the two-temperature source models is better con- 
strained by the SIS than by the LECS (see Fig. 5). Naively, 
one would have expected that the softer response of the LECS 
should have led to smaller uncertainty on the cooler tempera- 
ture. This example shows how the complex interactions of pa- 
rameters in the fit process can lead to non-intuitive behavior 
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and points, a posteriori, to the importance of basing the fit pro- 
cess on detailed simulation work. 

Simulations can only provide a hmited insight as they are 
done on the assumption that both the detector's caUbration 
and the emissivity model are not affected by significant errors. 
However, they also show that even relatively small amounts 
of systematic errors can have a strong influence in the de- 
rived best-fit parameters. The statistical error of a 10 000 counts 
coronal spectrum in the peak (i.e. around 1 keV) is ~ 2.5% per 
effective resolution element (obtained adding up the number of 
counts in the energy range corresponding to the FWHM reso- 
lution at 1 keV), for both the LEGS and the SIS. This level of 
statistical uncertainty is sufficient to induce a 50% spread in 
the best-fit metallicity for SIS data at the 90% confidence level. 
Thus, similar levels of systematic errors, if present in the re- 
gions in which the spectrum has a high signal-to-noise, could 
induce comparable systematic shifts in the best-fit metallicity. 

The analysis of real SIS spectra done in Sect. 6 points to- 
ward the likely presence of systematic effects due either to cal- 
ibration uncertainties or to problems in the emissivity models 
or to the simphstic model used to fit the data (or to a com- 
bination of all the above effects), as shown by the strong de- 
pendence of the best-fit metallicity on the energy range used 
in the fit. In particular, our results on the soft sources /3 Get 
and Gapella could indicate a lack of predicted model flux in the 
region around 1 keV. 

The present results point toward several items which should 
be considered critically when discussing coronal metallicities 
derived from low-resolution X-ray spectra. The most important 
are: 

- Galibration uncertainties are likely to induce systematic ef- 
fects in the determination of coronal abundances. As shown 
above, calibration uncertainties of a few percent could in- 
duce systematic shifts in the best-fit abundances of several 
tens of percent. Galibration uncertainties in the spectral re- 
gion in which the signal-to-noise of stellar spectra is higher 
(around 1 keV) are potentially the most detrimental. Future 
instruments, such as XMM and AXAF, should thus pay 
particular attention to caUbration in this spectral region (as 
well as in the region below 0.5 keV). 

- Uncertainties in the plasma emissivity codes are as criti- 
cal as detector calibration issues. Again, the region around 
1 keV is the one which should be the made as reliable as 
possible. Flux deficits in this region will have a very strong 
influence in the derived source parameters, not in terms of 
larger confidence regions, but rather in terms of systematic 
shifts of best-fit values, due to the pivoting role of the re- 
gion of the spectrum with the highest signal-to-noise ratio. 
This becomes even more critical if individual abundances 
are left free to vary, as elements like Ga and Ar, which have 
their K-complexes in the hard region of the spectrum, have 
their L-complex lines also in the ~ 1 keV region. Again, 
the fit process will be more influenced by the higher signal- 
to-noise L lines than by the better resolved but noisier K 
lines. Given the strong blending in the 1 keV region, uncer- 



tainties in few of the Unes have a ripple-down effect on all 
abundances determined in the spectrum. 
- Assuming too simplistic a model for the fit (as illustrated 
here by fitting power-law temperature distributions with 
two-temperature models) can also induce systematic shifts 
in the retrieved abundance. As the evidence from EUVE 
spectral analyses of coronal sources points toward complex 
plasma temperature distributions, it is not unlikely that two- 
temperature analyses of real coronal spectra are affected 
by similar effects. In this regard the LEGS wide energy 
range provides improved diagnostic capabilities, thanks to 
the presence of two pivoting regions, one which is line-rich 
and one which is almost line-free. 

A hint of the possible inadequacy of the spectral model 
being used is, in the case of the LEGS, the derivation of 
different best-fit abundances as a function of the adopted 
low-energy cut-off. As shown by the fit of the SIS data 
of Gapella and /3 Get discussed here, this may also be true 
with ASGA, although caUbration uncertainties in the softer 
detector channels may also be playing a role there. 
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